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ABSTRACT 

Massive black hole (BH) mergers can result in the merger remnant receiving a 
"kick," of order 200 km s _1 or more, which will cause the remnant to oscillate about 
the galaxy centre. Here we analyze the case where the BH oscillates through the galaxy 
centre perpendicular or parallel to the plane of the galaxy for a model galaxy consist- 
ing of an exponential disk, a Plummer model bulge, and an isothermal dark matter 
halo. For the perpendicular motion we find that there is a strong resonant forcing of 
the disk radial motion near but somewhat less than the "resonant radii" r R where the 
BH oscillation frequency is equal one-half, one-fourth, (1 /6, etc.) of the radial epicyclic 
frequency in the plane of the disk. Near the resonant radii there can be a strong en- 
hancement of the radial flow and disk density which can lead to shock formation. In 
turn the shock may trigger the formation of a ring of stars near r R . As an example, 
for a BH mass of 10 s M & and a kick velocity of 150 km s _1 , we find that the resonant 
radii lie between 0.2 and 1 kpc. For BH motion parallel to the plane of the galaxy 
we find that the BH leaves behind it a supersonic wake where star formation may be 
triggered. The shape of the wake is calculated as well as the slow-down time of the 
BH. The differential rotation of the disk stretches the wake into ring-like segments. 

Key words: galaxies: kinematics and dynamics — galaxiies: nuclei — galaxies: struc- 
ture 



1 INTRODUCTION 

Recent breakthroughs in numerical General Relativity have 
led to predictions of large recoil velocities of merged bi- 
nary supermassive black-holes (BH) as the binary radiates 
away linear momentum as gravitational waves during the fi- 
nal stages of merger (Gonzalez et al. 2007; Campanelli et al. 
2007; Lousto, Campanelli, & Zlochower 2009). Typical kick 
velocities are of order ~ 200 km s _1 (Bogdanovic, Reynolds, 
& Coleman 2007). One possible result of these mergers is 
to cause the resulting remnant black hole (BH) to oscillate 
through the stellar disk on timescales on the order of Gyr 
(Kornreich & Lovelace 2008, hereafter KL08; Blecha & Loeb 
2008; Fujita 2009). If the merged BH receives an impulse 
from the merger, some of that motion will be transmitted to 
the galaxy by dynamical friction. It is clearly of interest to 
know whether the motion of the BH through the disk gener- 
ates observable changes in morphology and dynamics in the 
galaxy. Binary black-holes have already been observed as a 
double nucleus in a quasar (e.g., Decarli et al. 2009). Fur- 
ther, Comerford et al. (2009) show that BH binaries resulting 
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from recent mergers are observable in the spectra of as many 
as 40% of Seyfert 2 galaxies. The possible observable effects 
of free and "wandering" BH merger remnants on disk galax- 
ies, however, have not yet been fully analyzed, de la Fuente 
Marcos & de la Fuente Marcos (2008) discuss the formation 
of stars in the wakes of runaway black holes ejected from the 
host galaxy. Here we consider more slowly moving BHs which 
remain bound to the host galaxy. 

This work first analyzes the case where the ejected BH os- 
cillates with angular frequency Q bhz through the galaxy cen- 
ter in a direction normal to the disk. Later, we discuss the case 
where the BH is ejected parallel to the plane of the galaxy. 
Sufficiently large gas accretion by the binary BH system is 
predicted to drive the orbital and BH spins into alignment 
normal to the plane of the galaxy (Bogdanovic et al. 2007). 
However, the amount is uncertain due to the uncertainty in 
the ratio of the two viscosity coefficients for the (r, z) and 
(r,4>) motion in the disk (e.g., Natarajan & Armitage 1999). 
This alignment favors BH ejection in the plane of the galaxy 
(Campanelli et al. 2007). The merging process preceding the 
ejection is assumed to be sufficiently slow that the host galaxy 
has returned to an axisymmetric equilibrium state. The BH os- 
cillations are observed in simulations to persist for many pe- 
riods (KL08; Blecha & Loeb 2008; Fujita 2009). Under these 
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conditions we find that there is a strong resonant forcing of 
the disk radial motion near radii tr where 2nClbhz = 0, r (rR), 
with n = 1, 2, .., where Q r is the radial epicyclic frequency. 
Near tr there can be a strong enhancement of the density 
in a ring of the disk and shock formation which may trig- 
ger the formation of a ring of stars. We note that models for 
the formation of observed "ring galaxies" assume the (single) 
passage of one galaxy nearly through the center of another 
(Theys & Spiegel 1976, 1977; Lynds & Toomre 1976). The 
gravitational interaction of the two galaxies causes the for- 
mation of a pronounced ring(s) of enhanced density (in one 
or both of the galaxies) where star formation is expected to 
be enhanced. 

For the case of BH ejection in the plane of the galaxy we 
find that the BH leaves behind it a supersonic wake where 
star formation may be triggered. The shape of this wake is 
calculated and the slow-down time of the BH is estimated. 
The differential rotation of the disk stretches the wake into 
ring-like segments. 

Section 2.1 describes the gravitational potential of the 
equilibrium disk galaxy at the time of the BH merger. This in- 
cludes an exponential disk of stars and gas, a Plummer model 
the bulge, and an isothermal dark matter halo. Section 2.2 
derives the dynamical equations of the system as driven by 
the radial force due to a BH oscillating perpendicular to the 
disk. The frequencies and amplitudes of the BH oscillations 
are taken from the iV-body simulations of KL08, Section 2.3 
discusses the forced oscillations of the gas disk, and §2.4 the 
possible parametric instability of the disk. Section 3 discusses 
the response of a galaxy disk to a BH ejected in the plane of 
the galaxy. Section 4 gives the conclusions of this work. 



2 THEORY 

2.1 Equilibrium 

The equilibrium galaxy is assumed to be axisymmetric and to 
consist of a thin disk of stars and gas and a spheroidal distri- 
butions consisting of a bulge component and a halo of dark 
matter. The model is the same as that used in the simulations 
of KH08. The model is similar to that of Fujita (2008, 2009). 
We use an inertial cylindrical (r, <j>, z) and Cartesian (x, y, z) 
coordinate systems with the disk and halo equatorial planes 
in the z — plane. The total gravitational potential is written 



$(r, 



(1) 



where $d is the potential due to the disk, <J>6 is due to the 
bulge, and $ h is that for the halo. The galaxy may have a 
central massive black hole of mass AI bh in which case a term 
<E>i,h = —GMbh/y/r 2 + z 2 is added to the right-hand side of 
(1). The particle orbits in the equilibrium disk are approxi- 
mately circular with angular rotation rate fi(r), where 



n 2 (r) = 



19$ 

r dr 



+ 



(2) 



The equilibrium disk velocity is v = rQ,{r)cj) . A central black 
hole is accounted for by adding the term Qjj h = GA'hh /r' A to 
the right-hand side of (2). 

The surface mass density of the (optical) disk is taken 
to be Ed = E<joexp(— r/rd) with E^o and r d constants and 
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Figure 1. Rotation frequency of the disk / = U/2n, one-half the 
radial epicyclic frequency f r /2 = Q r /4n, and one-fourth this value, 
/r/4, for the galaxy model described in §2.2. 



Md — 27rr-|Edo the total disk mass. The potential due to this 
disk matter is 

GM d 



*d(r,0) 



Td 



l R[h{R)Ki{R) - h(R)K (R)] 



and the corresponding angular velocity is 



l GM d 

9 »-3 



[Ia(R)K (R)- Ii(R)Ki (R)] 



(3) 



(4) 



where R = r/(2rd) and the I's and K's are the usual mod- 
ified Bessel functions (Freeman 1970; Binney & Tremaine 
1987, p. 77). We consider the mass of the disk stars and gas 
is M d = 2.8 x 10 10 M Q and r d = 3.5 kpc following KL08. For 
these values, v d = \fGMdfr~d = 186 km/s. 

The potential due to the bulge component is taken as a 
Plummer (1915) model 

GM b 



+ r 2 



2)1/2 



(5) 



where Mb is the mass of the bulge and r b is its characteris- 
tic radius (Binney & Tremaine 1987, p. 42). This component 
is similar to the spherical component considered by Fujita 
(2008, 2009). We have 

^2 GMb 



{r 2 b + r 2 ) 3 / 2 



(6) 



W e take M b = 10 10 M Q and r& = 1 kpc so that v b = 
y/GMb/n = 208 km/s again following KL08. 

The dark matter halo is assumed to have the isothermal 
distribution with the potential 



$ h = -v 2 h ln(r 2 c +r 2 + z 2 ) 



(7) 



where r c is the core radius of the halo, and Vh is the circular 
velocity at distances larger than r c . This potential is assumed 
to apply out to a large distance say 50 kpc. Equation (7) im- 
plies 



+ r 2 



(8) 



Representative values are vh = 250 km/s and r c = 2 kpc. The 
total dynamical mass of the galaxy model is Af to t ~ 6 x 10 11 . 
We do not attempt to model the dark matter distribution 
at small radii where ACDM simulations predict a cusp with 
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Figure 2. Vertical oscillation frequency f b hz = Q.bhz/'Zir as a function 
of the initial velocity of the BH from KL08 for M bh = 10 s M a . The 
line through the points is simply a smooth fit. 
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Figure 3. Dependence of the vertical amplitude of the BH motion 
z b ho on the initial velocity from KL08 for Mbh = 10 s Mq. The line 
through the points is simply a smooth fit. 
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(Navarro, Frenk, & White 1996). Observed rota- 



tion curves of spiral galaxies are however well fit by profiles 
Pdm oc (jb + r) 1 with tb is larger than 5 kpc (Salucci 2001; 
Salucci & Burkert 2000; Burkert 1995). 

Figure 1 shows an illustrative rotation curve fl(r) /2ty not 
including the contribution of the BH. The contribution of the 
BH for M bh = 10 s M is not significant for r > 0.2 kpc. 

2.2 Oscillating Black Hole 

We first consider the case where the BH of mass Mbh is ejected 
vertically and oscillates about the midplane of the galaxy, 



Zbh 



(9) 



where Q. b hz = 2-Kf bhz is the angular frequency of the vertical 
oscillation, and cs , c 5 , . . account for the fact that the motion 
is not in general simple harmonic, but the effective potential 
is an even function of z. However, for z m < 1 kpc, the mag- 
nitudes the c's is small compared with unity (KL08), and we 
neglect them. Figure 2 shows the dependence of f bhz on the 
initial kick velocity of the BH from KL08. We assume that the 
damping time-scale of the oscillations due to dynamical fric- 
tion is significantly longer than the oscillation period as found 
by KL08, Blecha & Loeb (2008), and Fujita (2009). 

The oscillating BH gives a time-dependent radial force or 
acceleration on the disk matter at a radius r, 



a r (r, t) 



GM bh r 



[r 2 + z! h (t))^ 



(10) 



where Q bh = y/GM bh /r 3 , and 

-3/2 



T = 



1 + 



Zbh 

r 



= go + gi cos(2Sl b hzt) + g 2 cos(4f2 Wiz t) + . 



(11) 



Here, go, g\, and g^ are dimensionless functions of z m /r as 
shown in Figure 4. The go term acts to modify the disk equi- 
librium mainly at distances <, 0.2 kpc, and it is not considered 
further. 

We model the disk as an axisymmetric thin fluid layer 
acted on by only the radial force of the BH (a r ). The disk 
consists of gas and stars so that its linear response to the 
black hole should be treated as two gravitationally coupled 




z m /r 



Figure 4. Dimensionless coefficients dependent only on z m /ro as de- 
fined in equation (11) and obtained numerically. 



fluids - gas and stars - with different surface densities (E 9 
and £*) and different 'sound' speeds (c s and c) following 
the approach of Lin & Shu (1970). In our galaxy E s <t; £* 
and c g <g c», The Toomre (1964) stability factor for axisym- 
metric perturbations of the gas Q g oc c 9 /E 9 is significantly 
less than that for the stars Q» oc c*/£„, A number of studies 
(e.g., Jog & Soloman 1984; Wang & Silk 1994; Rafikov 2001) 
) point out that the disk stability is determined mainly by the 
small-Q (less stable) component of the disk. For this reason 
we analyze the response of the gaseous component of the disk 
and let E 



du r 

~dt 
dt 



E 9 . We have 

vp ^4 

E r 



Or 



+ a r (r,t) , 



= 



(12) 



where d/dt = d/dt + u r (d/dr), E(r, t) is the surface mass- 
density and P(r, t) is the height-integrated pressure. The disk 
mass is conserved so that 



as 

dt 



1 d(ru r E) 
r dr 



= 



(13) 



There is also an equation for the gravitational potential. 

We linearize equations (12) and (13) by letting u r = 

+ 5u r , = rfl + Su^, P = P + SP, E = E + <5E, and 
$ = $0 + 5$. The disk equilibrium has = ~T, ~ 1 dP /dr + 
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Figure 5. Resonant radii where flr(rm) = 2Ubhz and U r (rR 2 ) 
iflbhz as a function of the initial BH velocity for M\,h = 10 8 Mq. 
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dt 
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Equations (12) then give 

<5E dP _ 1 dSP 
dr Eo dr 
d(r 2 Q) _ 



2m u . 



85$ 
Or 



+ a r 



dt 



d{r 2 n) 

dr 



(14) 



We assume P oc E r with T the adiabatic index. Thus 
5P = c 2 g 5T, with c g = (I\P/E) 1/2 the sound speed in the 
disk assumed independent of r. Also, we assume short radial 
wavelengths with \dSP/dr\ ^> \SP\/r and verify this later. 
The continuity equation gives d(ST,)/dt = —Ho(dSu r /dr). 
The equation for the gravitational potential for z — gives 
jfc r j<5<!> = — 27rG5E, where we have taken the perturbations 
to be proportional to exp(ik r r). 

We now multiply the equation for 5u T by (r/fi) and apply 
d/dt to it. Omitting the terms nonlinear in 5u r , this gives 



d 5u r 
dt 2 



r\2 f\ 

" + 2nisGZ T — + Q,' 2 
or 



da 



9Q r 2 

+ 



d 



ra r 



(15) 



dt r dr 

Here, f2 r = [r~ 3 d(r 4 0, 2 )/dr} 1/2 is the radial epicyclic fre- 
quency, s — sign(fe r ), and T < 1 is a "reduction factor" 
which is unity for |fc r zo| < 1 and decreases as |fc r zo| increases 
from unity where zo is the half-thickness of the disk (see Shu 
1984). 

With a r = 0, Q r =const, and 8u r oc exp(ik r r — wt) 
in equation (15), we recover the dispersion relation for ax- 
isymmetric perturbations of the disk, to 2 = Q 2 + (k r c g ) 2 — 
2nGT, T\k r \ (Toomre 1964). 

With a r given by equation (10) and included in equation 
(15), there are forced oscillations of the disk due to the da r /dt 
term at radii where tt r ~ 2Qthz (for the g± term in equation 
11), Q r m AQ bhz (for the g 2 term), etc. The radii where these 
resonances occur are shown in Figure 5. 

There can be a parametric instability of the disk due to 
the (Q./r)8u r d(ra r /Q,)/dr term at radii where Sl r £s Sl bhz 
(for the gi term), fi r w 2Qthz (for the g 2 term), etc. We first 
discuss the case of forced oscillations. 

2.3 Forced Oscillations of the Disk 

In the case of forced oscillations in equation (15) we have 
da r /dt = 2nQ. 2 h rg n Q,bhz sin(2nfif,h z t) for n = 1, 2, .. in equa- 



tion (11). At the resonant radius ru in the disk, fl r (rR) = 
2nQbhz- Note that Q. r varies slowly across the disk so that 
we can write D. 2 = £l 2 R [l — 2\n\{r — r^J/m], where 
Q. r R = Q r (rR) and r\ = (r/Q r )(dQ r /dr) which is found to 
negative for the considered conditions. We let Su r (r,t) = 
Su r (r) sin(2nQ,bhzt) so that equation (15) becomes 

d 2 
dx 2 



d 



la h bx 

dx 



Su r 



(16) 



where 



fio = fl(r R ) 



2ttGT, T 

CqQ.0 



r R Q, , K = — g n 



QrR 



^bh^rR 



Note that for the considered thin disks c g /uo <S 1 so that 
the a;— coordinate is an expanded version of r — r R . For 
the assumed conditions we find T = 1. We can write a = 
2(Q r R./Qo)(T/Q), where Q — Sl r HC 9 /(7rGEo) is Toomre's 
(1964) factor with Q > 1 disks being stable to axisymmet- 
ric perturbations. 

We can simplify equation (16) by letting 5u r /c s = 
U(x)F(x) and choosing F = exp(iaa;/2). Then, 

d 2 U 
dx 2 



a. 

— + bx 
4 



U — exp(— iax/2)K,(x) 



(17) 



A WKBJ solution of the homogenous part of the equation with 
U oc exp[i J x q(x')dx'\ gives q = (a 2 /4 + bx) 1 / 2 . Thus there 
is wave-like propagation for x > —a 2 /Ab which begins inside 
the resonant radius where x = 0. The region x < —a 2 /4b is 
forbidden and we assume U(x — > — oo) — > 0. 

The exact inhomogeneous solution to equation (16) is 



r„ iax/2 

OUr 7te ' 

c 



%y/2 



\A\{-a-px) / dyBi(-a - /3y)K(x)e 

■3 P 1 J xt 

-Bi(-a-/3a;) f dyAi(-a - Py))C{x)e- tay/2 } , (18) 

where a = a 2 /(4fe 2/3 ), /3 = b 1/3 , and Ai & Bi are the usual 
Airy functions. The value of xt is chosen to be well inside the 
forbidden region; that is, xe <§C — a 2 /4b. 

Figure 6 shows an illustrative solution of equation (16) 
for 5u r /c g and the corresponding fractional surface den- 
sity variations <5E/E. For this figure we consider a BH with 
M b h = 10 s Mq, a resonant radius t-r = 0.5 kpc where 
0. r (r R ) — 4Qbhz (that is, n = 2). This corresponds approx- 
imately to the results of KL08 with Vi n n w 150 km/s and 
z m ~ 1.3 kpc. From Figure 4 we find g 2 ~ 0.16. For the galaxy 
model of §2.1 we find fbhz = Slbhz/2-K w 16.9 Gyr -1 (or a 
period of 59.2 Myr), 27r/n rfi w 14.8 Myr, fi r fl/fi ~ 1-89, 
and rj w —0.311. We assume the galaxy disk of §2.1 has 
a gas mass-fraction of 0.22 and that c g /uo = 0.05 so that 
the Toomre factor of the gas disk, Q — 2.1 which corre- 
sponds to the disk being stable to axisymmetric perturba- 
tions. For the §2.1 model, uq = rRQ(rR) = 110 km s _1 so 
that c g = 5.5 km s _1 . We assume xe = —10. For these val- 
ues we find a = 1.8, b = 0.111, a = 3.52, /3 = 0.480, and 
K w -0.438[1 + x{cg/r R n )]- 3 . 

The full-width of the surface-density peak in Figure 
6 at half-maximum is Ax w 2. In terms of the ac- 
tual radius this translates to Ar « rR(c g /uo) 2 ^ 3 K , where 
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+ n^Ji - G n cos[2(a rR + 8n r )t]jsu r = o , (19) 



Figure 6. Illustrative solution of equation (16) for 5R(Su r )/c g (solid 
curve) and the corresponding 5S/So (dashed curve) for conditions 
described in the text. For this case x = 20(r — r R )/r R . Note that 
Su r (x, t) = 3t[6u r (x)] sin(U rR t) and 8T,(x,t) = 8T:(x) cos(U rR t) 
with <5E(ie)/£o = {£lo /£lrR)d\yt(&u r {x) / c g )}/ dx. The location of 
the peak of 8u r at x —5.4 coincides approximately with the loca- 
tion of the maximum of the Airy function Ai(— a — fix). 



K w {2n /Q. rR ) 2/ ' i \r 1 \- 1/3 « 1.53. We find Ar « 
0.1 kpc. This length is larger than the "Toomre-length" 
k^ 1 — r R (c g /uo)(Q,o/fl rR )Q, where kr is the "least stable 
wavenumber" (where uj(k r ) 2 is a minimum) for axisymmet- 
ric perturbations. For the conditions of Figure 6, fc^ 1 £s 0.028 
kpc. 

The magnitude of the disk response, 5u r /c g and <5E/Eo, 
to the oscillating BH is directly proportional to the BH mass 
Mbh and inversely proportional to the Q-factor of the gas disk. 
The response decreases strongly as the resonant radius r R in- 
creases, roughly as r^ 3 . 

A sufficiently strong disk response can lead to shock for- 
mation. The time-scale for a shock to form is of the order of 
t s h = \d6u r /dr\ . If t B h is less than half the period of os- 
cillation of 5u r , 7r/fi r , then a shock has time to form. Using 
the continuity equation, the condition for a shock to form is 
|5S|/So > 7t _1 . This condition is satisfied for a range of r 
for the conditions of Figure 6. The compression of gas in the 
shock wave may in turn lead the formation of a ring of stars. 

The amplitude of oscillation of the BH, z m , decays on 
a time-scale longer than its oscillation period (KL08) due to 
dynamical friction of the BH with the different components 
of the galaxy. Part of the decrease in the BH kinetic energy 
goes into the energy of the disk perturbation. The decrease of 
z m leads to an increase in the BH oscillation frequency flbhz 
(KL08). Thus the radius of a given resonance r R decreases 
with time. For the case of Figure 6, we estimate that shock 
formation occurs only for r R < 1 kpc, but continues down to 
r R ~ 0.2 kpc. 



2.4 Parametric Instability of the Disk 

Here we discuss briefly the influence of the term 
(Q/r)Su r d(ra r /0)/dr in equation (15) which can give rise 
to a parametric instability in the disk. This term can be writ- 
ten as G n n^5w r cas(2n£lbhzt), for n = 1, 2, .., where G n = 
—(Q,/r)Q~ 2 d(r 2 Q. 2 , h g n /Q)/dr is a dimensionless factor. The 
dominant parametric instability occurs when nfl b hz ~ £7r- 
Omitting for simplicity the r— derivative terms and the da r /dt 
term in equation (15) gives 



d 8u r 
dt 2 

where in this case the "resonant radius" r R is such that 
nftbhz = Qr(r R ), and 5Q r = rjQ, r Sr/r R is a measure of the 
radial distance from r R . The solution of this equation is oscil- 
latory with angular frequency Q, rR + 5Q r with an amplitude 
exponentially growing with growth rate 



G n 



V Qr R 



1/2 



(20) 



for JO r < Q rR \G n \/4, (Landau & Lifshitz 1960). The max- 
imum growth rate is max{uji)=Q rR \G„\/4. For the galaxy 
model of §2.1 and a BH mass M bh = 10 s M Q , we find that 
\G n \ decreases from roughly 0.15 at r rR — 0.2 kpc to 10~ 3 
at r r R = 1.5 kpc assuming g n — 0, 16. We conclude that 
the parametric instability is unimportant compared with the 
forced motion (§2.3) for r rR > 0.2 kpc. 



3 BH EJECTION PARALLEL TO THE PLANE OF THE 
GALAXY 

For the case where the BH is eject parallel to the plane of the 
galaxy we can assume that the motion is in the x— direction, 



Xbh = Xm Sm(Qbhrt) , 



(21) 



with x m the amplitude of the motion and ftbhr the angu- 
lar frequency as determined by KL08. The motion of the BH 
through the gas disk of the galaxy is highly supersonic. Thus 
the BH will leave behind it a narrow conical shock wave 
or wake. This shock wave may trigger star formation par- 
ticularly in regions of high density such as dense molecular 
clouds. When the BH is at say Xbh(t') it leaves a "shock- 
section" at the radius |a;6/i(t')| . At a later time t > t' this 
shock-section rotates about the galaxy's centre by an amount 
A(f> = Q,(\xbh\){t — t') > 0. following the rotation curve of the 
galaxy. Figure 7 shows the geometry of the wake for the case 
where the BH has gone through one period of radial oscilla- 
tion, 2-rr/Qbhr- For the case shown this period is 59 Myr. The 
differential rotation of the disk stretches the wake into ring- 
like segments of radii approximately equal to the maximum 
excursion of the BH, x m , where the BH most slowly. 

The relative velocity of the BH and galactic gas disk is 
Av = <p u^(r) — xdxbh/dt , where u^s is the rotation velocity 
of the disk. As mentioned the flow is highly supersonic with 
|Av| S> c g , where c 9 is the sound speed in the gas. In this 
limit there is Bondi, Hoyle, & Lyttleton accretion to the mov- 
ing BH (Hoyle & Lyttleton 1939; Bondi & Hoyle 1944) where 
gas approaching the BH with an impact parameter less than 



ba = 



2GM bh 
|Av| 2 



(22) 



accretes to the BH if b a < h with h the half-thickness of the 
disk. In this case the cross-section for accretion is ivb 2 . If b a > 
h the cross-section is approximately Ab a h. We assume b a < 
| Xbh | during most of the BH's orbit. The relevant case is b a > h 
and this gives the an accretion rate M b h = 4:b a hp g \ Av| where 
p g is the gas density. The accretion time-scale is then 



T — 



Mb, 
Mb, 



|Av| 
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Figure 7. Wake of the BH ejected in the plane of the galaxy for the 
case where the BH has gone through one period of radial oscillation, 
27r/f!6h r . For this plot x m = 1 kpc and f^ r = = 17 

Gyr -1 corresponding to a period of 59 Myr. The points marked by 
a, b, .. correspond to the t' = 0, t' = ir / fitter) , etc. 
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For these reference values b a « 87 pc and M bh « 
O.43M yr _1 . The drag force on the BH is simply Fth = 
ip g b a h\ Av| Av. This acts to both slow down the BH radial 
motion and impart to it azimuthal motion in the rotation di- 
rection of the galaxy. The time-scales for the changes in the 
motion are all of order T a ■ 



4 CONCLUSIONS 

Here, we first analyzed the linear axisymmetric perturbations 
of a gas disk driven by a black hole oscillating vertically along 
the axis of symmetry. We find that there is a strong reso- 
nant forcing of the disk radial motion near "resonant radii" 
tr where the BH oscillation frequency is equal one-half, one- 
fourth, (1/6, etc.) of the radial epicyclic frequency in the 
plane of the disk. Near the resonant radii there can be a 
strong enhancement of the radial flow velocity and disk den- 
sity which can lead to shock formation. The shock formation 
occurs during one period of oscillation of the BH which is 
assumed longer than the BH damping time due to dynami- 
cal friction. This shock may trigger the formation of a ring 
of stars near r R . As an example, for a BH mass of 10 8 M 
and a kick velocity of 150 km s _1 , we find that the resonant 
radii lie between 0.2 and 1 kpc. The magnitude of the disk re- 
sponse is proportional to the BH mass, inversely proportional 
to the Toomre Q of the gas disk, and decreases rapidly as m 
increases. The resonant radii increase as the initial BH kick 
velocity increases (KL08). 

For BH motion parallel to the plane of the galaxy we find 
that the BH leaves behind it a supersonic wake which over 
time gets contorted into a complicated shape by the galaxy's 
differential rotation. The shape of the wake is calculated for 
an illustrative case as well as the slow-down time of the BH. 



The differential rotation of the disk stretches the wake into 
ring-like segments of radii approximately equal to the maxi- 
mum excursion of the BH, x m , where the BH most slowly. 

Many other processes may be involved in the formation 
and evolution of nuclear rings in galaxies (e.g., van de Ven & 
Chang 2009). 
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